Variations on the adiabatic invariance: the Lorentz pendulum 
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We analyze a very simple variant of the Lorentz pendulum, in which the length is varied exponentially, 
instead of uniformly, as it is assumed in the standard case. We establish quantitative criteria for the condition of 
adiabatic changes in both pendula and put in evidence their substantially different physical behavior with regard 
to adiabatic invariance. 



I. INTRODUCTION 

As early as 1902, Lord Rayleigh 1 investigated a pendu- 
lum the length of which was being altered uniformly, but 
very slowly (some mathematical aspects of the problem were 
treated in 1895 by Le CornvP). He showed that if E(t) denotes 
the energy and £(t) the length [or, equivalently, the frequency 
v(f)] at a time f, then 
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E(0) 
v(0) 



(1.1) 



This expression is probably the first explicit example of an 
adiabatic invariant; i.e., a conservation law that only holds 
when the parameters of the system are varied very slowly. The 
name was coined by analogy with thermodynamics, where 
adiabatic processes are those that occur sufficiently gently. 

At the first Solvay Conference in 191 1, Lorentz, unaware of 
Rayleigh's previous work, raised the question of the behavior 
of a "quantum pendulum" the length of which is gradually al- 
tered^ (by historical vicissitudes 4 his name has become inex- 
tricably linked to that system). Einstein's reply was that "if the 
length of the pendulum is changed infinitely slowly, its energy 
remains equal to hv if it is originally hv", although no detail 
of his analysis are given. In the same discussion, Warburg in- 
sisted that the length of the pendulum must be altered slowly, 
but not systematically. As Arnold aptly remarks^ "the person 
changing the parameters of the system must no see what state 
the system is in. Giving this definition a rigorous mathemat- 
ical meaning is a very delicate and as yet unsolved problem. 
Fortunately, we can get along with a surrogate. The assump- 
tion of ignorance of the internal state of the system on the part 
of the person controlling the parameter may be replaced by 
the requirement that the change of parameter must be smooth; 
i.e., twice continuously differentiable". This important point 
is ignored in most expositions of adiabatic invariance. 

Ehrenfest, who did not attend the Solvay Conference 
and was not cognizant of that discussion, had indeed read 
Rayleigh's paper and employed those ideas to enunciate his 
famous adiabatic principle 6 , which was promptly reformu- 
lated by Born and Fock 7 in the form we now call adiabatic 
theorem 8 . In fact, this was a topic of uttermost importance 
in the old quantum theory 9 '°. To put it simply, if a physical 
quantity is going to make "all or nothing quantum jumps", it 
should make no jump at all if the system is perturbed gently, 
and therefore any quantized quantity should be an adiabatic 
invariant. The reader is referred to the book of Jammer-^ for 
a masterful review of these questions, as well as the lucid and 



concise mathematical viewpoint of Levi- CivittaEEl. 

The topic of adiabatic invariance has undergone a resur- 
gence of interest from various different fields, such as plasma 
physics, thermonuclear research or geoph' L sicP, although 
perhaps Berry's work on geometric phases^ has put it again 
in the spotlight. The two monographs by Sagdeev et ai 16 and 
Lochak and Meuniei 17 reflect this revival. More recently, the 
issue has renewed its importance in the context of quantum 
control (for example, concerning adiabatic passage between 
atomic energy level a 18 " 23 ! ), as well as adiabatic quantum com- 
putatiorP4n2D. 

Adiabatic invariants are presented in most textbooks in 
terms of action-angle variables, which involves a significant 
level of sophistication. Even if a number of pedagogical pa- 
pers has tried to alleviate these difficultie J 27 " 34 !, students often 
understand this notion only at a superficial level. Actually, 
perfunctory application of the adiabaticity condition may lead 
to controversial conclusions, even in the hands of experienced 
practitioners^N22l. 

Quite often the Lorentz pendulum is taken as a t ypica l ex- 
ample to bring up this twist for graduate students^E2l j n 
spite of its apparent simplicity, the proof of invariance is gen- 
uinely difficult 41 48 and details are omitted. The purpose of 
this paper is to re-elaborate on this topic, putting forth per- 
tinent physical discussion that emphasizes the motivation for 
doing what is done, as well as to present some variation of the 
Lorentz pendulum in the hope that its solution will shed light 
on the subject at an intermediate level. 



II. THE UNIFORMLY VARYING PENDULUM 

A. Basic equations of motion 

We confine our attention to the ideal case of a simple pen- 
dulum of mass m and variable length £(t), oscillating under 
the gravity. The Lagrangian of the system is 



1 

L = —m 
2 



dt 



-mgicosd- , (2.1) 



where #(?) denotes the inclination of the pendulum with the 
vertical. The Euler-Lagrange equation for the generalized co- 
ordinate # becomes 



d 2 $ 
dt 2 



dV 
dt 



sin# =0. 



(2.2) 
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Note in passing that the length £(t) acts as a geometrical 
(or holonomic) constraint, which here becomes time depen- 
denPSD. 

In what follows we shall restrict ourselves to the regime 
of small oscillations (that is, sin # ~ In this Section, we 
deal with the example of a pendulum for which the length is 
uniformly altered in time; i. e., 

£(t) = £ {l + et), (2.3) 

where e is a small parameter with the dimensions of a recip- 
rocal time. It will be convenient to let T = et and define a 
dimensionless time-dependent frequency 
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Equation (2.2 1 can be thus recast as 



0- 



-2-# + cd 2 (t)#=0, 



(2.4) 



(2.5) 



where the dot represents differentiation with respect to T. We 
limit our analysis to a lengthening pendulum because if the 
amplitude of the initial displacement is small, then the result- 
ing displacement will stay small. On the other hand, if the 
pendulum is shortening, then even if the initial displacement 
is small, the subsequent displacement will grow in time, vio- 
lating the linearization hypothesis. 

Using basic properties of Bessel functions^, the following 

= -7= = [A Ji (2(HdVT+t) +B Yi (2cjdVT+t)] , 
Vl +t 

(2.6) 

is a solution to (2.5), with (Oq = ffl(O) and J n (x) and Y„(x) 
the Bessel functions of nth order and first and second kind, 
respectively. The constants A and B must be determined by the 
initial conditions, which we take, without loss of generality, as 
t?(0) = #0 and t?(0) = 0. The final result reads as 



d(T) = 



[/ 2 (2cao) ri(2fflbvT+T) 



VT 

y2(2a>o)./i(2aWl + T) 



(2.7) 

As a side comment, we remark that, in spite of its usual des- 
ignation, neither Rayleigh nor Lorentz actually examined the 
behavior of this pendulum; this was first ac comp lished much 
later by Rrutkov and Fock 41 , who obtained (2.7 1 and also de- 



rived Eq. (JTTTJ directly therefrom. Indeed, this solution allows 
one to investigate in great detail the periods of this systerrPO. 



Since e appears in the denominator in the definition (2.4 1, 
<X>o is actually very large (for example, if £q = 1 m and e = 
0.01 s , then C0q ~ 313). This suggests to consider the limit 



2coo\/l + T> 1 . 



(2.8) 



and then take the leading term in the asymptotic expansion of 
the Bessel functions^!] 



J n (x) ~ \l — COS \x- - - — j , 



(2.9) 
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FIG. 1. Exact solution (red continuous line) and asymptotic ap- 
proximation (blue points) for the uniformly varying pendulum with 
e = 0.1 s" 1 , t = 1 m and i9q = 0.3 rad. 



Consequently, we get 
J?0 



d(T) 



(l+T) 3 /4 



COs[2(Bt)(v / l +T- 1)] ■ 



(2.10) 



By making use of the approximation (2.10 1, one can check 
that the maximum angular amplitude # max scales as 



$max(^) 
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(2.11) 



which shows that it is a decreasing function of time. More- 
over, lim T ^ OT #max(T) = 0. 



In Fig. 1 both the exact solution (2.7i and its asymptotic 



approximation (2. 10 1 are plotted. The error associated with 
( 2.1 0| > is very small; obviously, this error is smaller when T 
is larger, since larger t's improve the approximation of ( |2.8| l. 
This can be formally expressed as 



lim 2coo\/l + T = 00 • 



(2.12) 



This limit guarantees that, for a given value of Oo, the asymp- 
totic expansion is always a good approximation for any T. 
From this perspective we can assert that the validity of the 
asymptotic expansion depends on the values of #0 (and e) but 
it is independent of time. 



B. Adiabatic invariance 

As pointed out in the Introduction, the concept of adiabatic 
change is associated with a variation that occurs infinitely 
slowly. The observer who is controlling the changes does not 
know the internal state of the system. In practical terms, this 
means that the change is adiabatic when the variation is car- 
ried out continuously and so slowly that the change 8£ of the 
length is very small compared to the length £ of the pendu- 
lum^, i.e. 



St 



< 1. 



(2.13) 
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By considering that the temporal interval in which the varia- 
tion 81 is produced coincides with the local period T for the 
length £, we have 



8£: 



dt 1 



(2.14) 



and recalling that T 
example at hand as 

U 



2%^J £jg, we can rewrite (2.13 1 for the 0.86 
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(2.15) 



with Tq being the period of the pendulum at T = 0. This re- 
quirement is independent of j9q and, thus, independent of the 
amplitude of the oscillations, which is intuitively expected. 



Equation ( 2. 15 1 clearly suggests that if the change in length 



is initially adiabatic, it will remain forever. Moreover, the adi- 
abatic character of the system will improve as time goes on. 



Interestingly enough, condition (2.15 1 is formally equiva 



lent to (2.8 1 ensuring the validity of the asymptotic approxi- 
mation; they will always be satisfied whenever 



(2.16) 



which, from the definition of e, can be equivalently recast as 



£ <S Clim = 



(2.17) 



Equations ( 2. 16 1 or ( 2. 17 1 (which do not depend on time) pro- 
vide a sensible criterion for the adiabatic change in a uni- 
formly varying pendulum. The lesser the initial length, the 
more quickly can be lengthening the pendulum under the adi- 
abatic hypothesis. Alternatively, <»ij m can be seen as the mini- 
mum value of ©o for which the asymptotic expansion is valid 
independently of T. 

To give a more quantitative argument, we compute the func- 
tion 



7(T) = 



V(T)' 



(2.18) 



that turns out to be an adiabatic invariant for arbitrary periodic 
motions in one degree of freedorrP. Here 77(f) is the Hamil- 
tonian and v(t) the frequency of the oscillations. 
For a pendulum, the Hamiltonian is 



77 



2 



gW 2 



(2.19) 



where p = £ 2 d~&/dt is the generalized momentum conjugate 
to r>. 

Taking into account the relations^ 

^-[x- l J x {x)] =x- 1 J 2 (x), -rix^Yiix)] =x- 1 Y 2 (x), 
ax ax 

(2.20) 

and the solution ( |2.7[ ), the angular velocity of the pendulum 
can be expressed as 

i? = ,M2(2cooVT+^) , (2.21) 




FIG. 2. Plot of 7(t) as a function of T for the uniformly vary- 
ing pendulum with £q = 1 m and $o = 0.3 rad. The curves corre- 
sponds to (Oq — 2, 10 and 1000 (which are associated with the values 



e = 0.2491,0.0498 and 0.00005 s" 
plitudes. 



), following the decreasing am- 



which immediately leads to 

H{%) = H Q n 2 (ol[J^ 2 \{2ox ) y/l + x) + Jt% (2&W1 + t)] . 

(2.22) 

Here, for notational simplicity, we have introduced the func- 
tions 

Jfti (x) = / 2 (2fflo) Yi (x) - Y 2 (2coo) /, (jc) , 

(2.23) 

#fa(x) = Y 2 (2coo) J 2 (x)~J 2 (2coo) Y 2 (x) , 

and Hq = $q(Dq£ 2 £1/2 is the total energy at T = 0. Finally, 
since v(t) = Vo/y/1 + T, with vo = (Qo£/2n, we get 

^0 f^^T„2 



7(t) 



Vo 



x [Jf 2 l(2(0 O VT+T) + J^ 2 ]{2(OoVT+T)}. (2.24) 

Because 7(t) is a time-dependent function, it will not be, in 
general, an adiabatic invariant. In other words, for arbitrary 



values of the parameters, the solution (2.7i will not be as 



sociated with adiabatic changes in the length of the pendu- 
lum. The function 7(t) is shown in Fig. 2 for different val- 
ues of COq. One immediately concludes that the larger Oo, the 
lesser time-dependent 7(f) becomes, which is in full agree- 
ment with ( |2.16| l. 

To complete the analysis, we proceed to calculate 7(t) with 
the asymptotic approximations for the Bessel functions as in 



Eq. (2. 10 1. After a direct manipulation, we end up with 



77(T) 



77 



so that 



Vo 



(2.25) 



(2.26) 



This is an important result: for large enough values of t, the 
quantity 7(t) becomes an adiabatic invariant. This validates 



4 



the previously suggested conclusion: the condition establish- 
ing the validity of the asymptotic expansion of the Bessel 
functions is conceptually equivalent to the condition of adi- 
abatic invariance. 



III. THE EXPONENTIALLY VARYING PENDULUM 

A. Basic equations of motion 

We turn now our attention to the instance where the length 
of the pendulum is altered non uniformly. More concretely, 
we take 



l(t)=tot 



(3.1) 



Equation (2.2 1, when small oscillations are considered, re- 
duces in this case to 



ti+2$ + coie- z $ = 0. 
The change # = 0e~ T gives 

d + {a>le- x -l)d = 0, 



(3.2) 



(3.3) 



which has again an exact solution in terms of Bessel func- 
tions^. The result, employing the original variables, reads as 



0(t) = e~ r [A J 2 (2cOoe 



-t/2\ 



BY 2 (2cOoe- T / 2 )}. (3.4) 



To fix the constants A and B we take the same initial conditions 
as before, namely #(0) = t9o and #(0) = 0. Applying the 
relations 

~[x 2 J 2 (x)]=x 2 J 1 (x), ~[x 2 Y 2 (x)]=x 2 Y 1 (x) (3.5) 
ax ax 

in conjunction with the change of variable x — 2<»o<?~ T//2 and 
the Wronskian 



J 1 (x)Y 2 {x)-J 2 (x)Y 1 (x) = - 



Tlx 



(3.6) 



one can show that the final solution is 

0(t) - n$ ( m )e - z [Fi(2fflo) h{2(o a e- T l 2 ) 
-h{2(Oo)Y 2 {2(Doe- T l 2 )]. 



(3.7) 



Much in the same way as for the pendulum with uniformly 
varying length, if 



2(Doe 



-T/2 



> 1. 



(3.8) 



is satisfied, we can replace equation (3.7i by its asymptotic 
approximation, leading to 



0(t) ~ tV~ 3T/4 cos[2o)o(l -e~ T/2 )] 



(3.9) 



Within this approximation the maximum angular ampli- 
tude # max scales exactly as in Eq. ( |2.11| i, and also 
lim T ^ m # m ax( T ) = 0. However, as expected, the angular am- 
plitude # max falls off more quickly in this case. This behavior 




FIG. 3. Plot of the exact solution (in red) and the asymptotic approx- 
imation (blue points) for the exponentially varying pendulum with 
the same parameters as in Fig. 1 . 



can be clearly appreciated in Fig. 3, where both the exact solu- 



tion (3.7 1 and its asymptotic approximation (3.9 1 are plotted. 
The agreement is again remarkable. 

In contradistinction with the situation described by 



Eq. (2.12 1, the exponentially varying pendulum leads to the 



limit relation 



lim 2a>oe 



-T/2_ 



= 0. 



(3.10) 



This points out the more important conceptual difference be- 
tween these two cases: for the uniformly varying pendulum 
the validity of the asymptotic expansion of the Bessel func- 
tions only depends on the values of e and #o, but it is inde- 
pendent on the time. For the exponentially varying pendulum 
the validity of that approximation is time dependent and for 
large values of the time this approximation breaks down. 

At first glance, this important difference between the two 
penduli can seems only a formal one. However, as we shall 
see in the following it has strong implications for the adiabatic 
invariance. 



B. Adiabatic invariance 

We next analyze the adiabatic change for this example. 
Equation (2.14i applied to ( |3.1| l gives as a requirement for 
adiabatic invariance 



81 



eT e 



T/2 



27t 

(Oq 1 



*/2 



«1, 



(3.11) 



which again is formally equivalent to the condition (3.8 i for 
the validity of the asymptotic approximation of the Bessel 
functions. 

This indicates that even if the change of length is initially 
adiabatic, it will remain so only for a finite interval of time. 
To put in another way, ( |3.1 1\ holds true whenever 



COq > e zl2 0) {{m - 



(3.12) 
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FIG. 4. Plot of /(t) as a function of X for the exponentially vary- 
ing pendulum with £q = 1 m and t?q = 0.3 rad. The curves corre- 
sponds to Wq — 2, 10 and 1000 (which are associated with the values 
e = 0.2491,0.0498 and 0.00005 s" 1 ), following the decreasing am- 
plitudes. 



The function 7(t) is represented in Fig. 4. As we can see, 
the fluctuations of 7(t) increase with time. Thus, for large 
enough time, 7(t) will never be an invariant quantity, irrespec- 
tive of the value of (Oq. However, for the time window chosen 
in the figure, we see that for C0o = 1000, 7(t) looks invariant 
over the entire interval. 

Our last step is to calculate 7(t) using the asymptotic ap- 
proximations for the Bessel functions. Now, we have 

H(T)^H e-^ 2 (3.17) 

so that 

/(T)^=,^ 2 tf o y /2 , (3.18) 

which is identical to what we have obtained for the uniformly 
varying pendulum. 

IV. CONCLUDING REMARKS 



or, recalling the definition of £»(t), 



21n 



eT 



(3.13) 



Accordingly, there does not exist a minimum fixed value of 
©o such that adiabaticity holds true forever. In our opinion 
this is the more important lesson from this paper: the adiabatic 
condition does not need to hold, in general, for all times. 
The formal equivalence discussed so far provides a math- 



ematical interpretation fo (3.13 »: the change will be adia 



batic in those conditions in which the asymptotic expansion 
of Bessel functions is justified. 

Note that the change of the length in this example is in- 
finitely continuously differentiable, but the adiabaticity holds 
true only in a time interval fixed by the values of £o and e. 
This by no means contradict Arnold's adiabaticity require- 
ment mentioned in the Introduction (that is, that the change 
of the length must be twice continuously differentiable), since 
Arnold explicitly considers a finite time interval in the defini- 
tion of the adiabatic invariants^. 

Finally, we calculate explicitly the total energy for this case. 



Using again the Hamiltonian (2. 19 1 and the solution (3.7 i and 
its time derivative, we get 



H{x)=Hq n 2 (0^e- z 

x [^ 2 1 (2ffl e- T / 2 ) 



^f 2 (2(Oo e - T/2 )], (3.14) 



with 



Jf u (x) = Fi(2ffl ) J 2 (x) - Ji(2fflb) Y 2 {x) . 



(3.15) 



MH i (x) = J x (2fflb) Yi (x) - Yy (2fflo) Ji (jc) . 
In consequence, 7(t) becomes 

/( T ) = — ^ 2 CO V T / 2 [^ 2 (2co 0e - T / 2 ) + ^ 2 2 (2to e- T/2 )]. 
vo 

(3.16) 



We have explored in detail two nontrivial yet solvable ex- 
amples of penduli of varying length with the purpose of a bet- 
ter understanding of the concept of adiabatic invariance. 

The ambiguous criteria of "infinitely slow variation", or the 
assumption of "ignorance of the internal state of the system on 
the part of the person controlling the variable parameter", usu- 
ally employed to establish the condition of adiabatic change, 
are replaced here by more quantitative criteria. 

For the two penduli considered in this paper, we have shown 
that the physical meaning of adiabatic change is formally con- 
tained in the mathematical condition of validity for the asymp- 
totic expansion of the Bessel functions: the validity of the 
asymptotic approximation implies adiabatic change and vice 
versa. 

The analysis carried out in this work invites to a more gen- 
eral reflection: it is important to pay special attention to the 
meaning of mathematical approximations. Actually, the ra- 
dius of convergence of some systematic approximation to an 
exact solution has always a physical origin. 
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